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By adding a linear term to a renormalization-group equation in a system exhibiting infinite-order 
phase transitions, asymptotic behavior of running coupling constants is derived in an algebraic 
manner. A benefit of this method is presented explicitly using several examples. 

PACS numbers: 64.60.Ak, 05.70.Fh, 05.70.Jk, ll.10.Hi, 02.30.Hq 



I. INTRODUCTION 

Renormalization-group (RG) technique is one of the most powerful methods for investigating critical phenomena in 
statistical physics Q. In general, RG transformation (RGT) consists of a coarse graining and a rescaling. It reduces 
many-body effects in a statistical model to an ordinary differential equation of coupling constants. The differential 
equation is called the RG equation (RGE), and has generally the following form: 



dg 

dt 



= V(g), 



(1) 



where g = (<?i,. .. ,g n ), & collection of coupling constants depending on t, and t = logL with L giving the length 
scale of the coarse graining in the RG. One obtains a beta function V (g) = (Vi(g), . . . , V n {g)) by applying the 
RGT explicitly to a statistical model. We can derive universal exponents that characterize critical phenomena from 
asymptotic behavior of solutions of Eq. |TIJ for large t. 

Since the asymptotic behavior is determined by vicinity of a fixed point g* , linearization of V (g) about g* is 
effective enough to obtain the exponents. For example, in a second-order phase transition, the correlation length £ 
typically behaves as 



£ = const./|T-T c 



(2) 



where v is the correlation-length exponent and T is a parameter specifying a state in a statistical model (e.g., the 
temperature). In the language of RG, T parametrizes initial values of RGE. The trajectory starting from the initial 
value at T — T c is absorbed into the fixed point. Other trajectories approach g* once but leave the fixed point 
subsequently, as shown in Fig. This implies that it takes longer for g to leave the fixed point as T approaches T c . 




FIG. 1: Typical RG trajectories near a phase transition. As T changes, an initial value moves on the dashed line. The trajectory 
with T = T c is absorbed into the fixed point. 
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If the scaling matrix M (g*), where 

Mij (g*) = ^(g*) (3) 

has a unique positive eigenvalue a, then this period behaves as const, x (T — T c )~ a . The universal exponent v is 
obtained from a by v — 1/a. Thus we do not need to find an explicit solution of generally nonlinear RGE JT]l. 
On the other hand, in the case of infinite-order phase transitions, £ behaves as 

£ = const, x exp (A/\T - T c \ a ) , (4) 

where a is a universal exponent and A is a non-universal constant. Such behavior is observed when all the coupling 
constants are marginal, i.e., the canonical dimensions of the coupling constants are zero at g*. Since the linear term 
in V (g) is proportional to the canonical dimensions of g, Mij(g*) = for all i and j. It indicates that we cannot 
extract the asymptotic behavior from the scaling matrix M(g*) in an infinite-order phase transition, in contrast to 
a second-order one. Therefore, explicit solutions were traditionally required in the case of an infinite-order phase 
transition such as the BKT phase transition^- 

This difficulty has been recently overcome in Ref.0, where an RG for RGE JTJ is used for deriving asymptotic 
behavior of solutions. A general idea of RG, applied as a tool for asymptotic analysis of non-linear differential 
equations, is developed in Refs.j3,0|- 

In this report, we present another method. Namely, we derive a from the following deformed RGE: 

^ = e(g-g*)+V(g) = V(g), (5) 

where e is a real number but not necessarily small. As we will see in the next section, the RG equation (JJJ for the 
RGE JIJ has a complicated form compared with the deformed RGE. Hence, using the deformed RGE makes derivation 
of the critical exponent simple. Another benefit of this approach is as follows: suppose that an infinite-order phase 
transition occurs when the spatial dimensions of the original statistical model are d c . Then, the deformed RGE can be 
derived when they are d c — e, under the condition that all the coupling constants have a common canonical dimension. 
This condition is satisfied by various field-theoretical models, e.g., an effective theory of antiferromagnets0, a model 
containing several gauge fields0, a model describing true self-avoiding random walksQ, and a model of nematic 
elastomers In Ref.Q, infrared asymptotic behavior in d c dimensions and that in d c — e dimensions are analyzed 
separately because of the problem of the vanishing scaling matrix explained above. Our method enables us to obtain 
universal quantities in both of the cases simultaneously. We will show this advantage in the last example of Sec II VI 



II. RGE FOR RGE 

Here we summarize the results of Ref.|3( that will be used later. We consider an RGE (JIJ for infinite-order phase 
transitions that are controlled by a fixed point g* . In what follows, we put g* = for convenience. Suppose that we 
have obtained V(g) by the lowest-order perturbation. Since linear terms vanish in infinite-order phase transitions, 
components of V(g) are quadratic in g. Hence the scaling property 

V(kg) = k 2 V(g) (6) 

holds in this case. The algebraic method to compute a shown in Rcf. 3] employs another RG to Eq. which is 
defined as follows: let g(t, ao) be the solution of Eq. (JTJ with the initial condition g(0, ao) — Uq. Choose a real number 
r and we evolve g in time by s(r), such that e T g(s(T),ao) £ 5*, where S is a sphere with the radius ao = |ao| and 
with the center at the origin. Thus we have the map R T : ao i— > e r g(s(r), ao) = a (r). Thanks to Eq. ©, R T satisfies 
the semigroup property R T + T ' = R t Rt', so that R T is called RG for RGE (JIJ. The infinitesimal transformation leads 
to the following new RGE: 

da a( „s- p ( a ) V ( a ) 2 ,„s 

— =P{a) = — a , 7 

dr a ■ V (a) 

where P (a) is the projection operator defined by Pij (a) = Sij — djOj/oQ. Since (3 (a) is perpendicular to a, solutions 
of the new RGE are restricted on S. Introducing the polar coordinates {# a }i<a<n-i on S and the corresponding 
orthonormal basis, 



e a = f a (a) 1 f a (a) = 



da 



(8) 
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FIG. 2: Schematic trajectories of RGE. The solid lines are for the original RGE Q while the dashed lines are for the new 
RGE Q defined on S. Here (a) is the case where a unique positive eigenvalue exists in n(a*). (b) is the case where all the 
eigenvalues of /J.(a*) are negative. 

we can expand /3(a) as 

n-1 



/3(a) = Y,P a (a)e a . (9) 



a=l 



The new RGE is represented as 



-^{a) = f- 1 {a)~p a {a) (10) 



by the polar coordinates. 

It is easily found that a* G S is a fixed point of the new RGE (JJJ if g(t, a*) is a straight flow line. In particular, 
a fixed point on an incoming straight flow line satisfying a* ■ V(a*) < plays an important role because trajectories 
near this fixed point correspond to trajectories of Eq. approaching g*. In contrast to the original RGE, we can 
generally linearize the new RGE about a*. In Rcf. 3], it is shown that the scaling matrix of the new RGE 

MQ/3 (a*)^/- 1 (a*)|^(a*) (11) 

plays a similar role to M(g*) in the original RGE describing a second-order phase transition. Namely, if the matrix 
/i(a*) has a unique positive eigenvalue A, in which typical trajectories of the original RGE are in Fig. [21 (a), we can 
observe divergence of the correlation length by one-parameter tuning and 

' " 5 < 12) 

in Eq. J2J. On the other hand, if all the eigenvalues of \i (a*) are negative, where typical trajectories are in Fig. [21 
(b), g(t,a ) behaves as 

9{t > ao) ~cWj-t e *- (13) 

In this formula, e* = a*/ao and C (g) is defined by the relation 

C(g)\g\ 3 = -g-V(g). (14) 

The asymptotic behavior in Eq. I|13(l is important for investigating finite size scaling in a statistical system, for 
example. 



III. DEFORMED RGE 



Next, we consider the deformed RGE (0 putting g* = 0. We can take e > without loss of generality. A fixed 
point c* of the deformed RGE solves 



V(c*) = ec* + V (c*) = 0. 



(15) 



(a) (b) 
FIG. 3: (a) Schematic trajectories for the original RGE. (b) Those for the deformed RGE. 



A key feature of the deformed RGE is that c* in Eq. Ijl5|l and a fixed point a* of the new RGE JJJ) on an incoming 
straight flow line has one-to-one correspondence via 

a* = ^c* (16) 

c* 

as depicted in Fig. |3] Writing V (g) as 

71-1 

V(g) = J2V<x{g)e a + V n (g)e n , (17) 

Q = l 

where e„ = g/g, we have the deformed RGE in the polar coordinates: 

^-(g) = f* 1 (g)V a (g) 

J (9) = eg + V n (g) . (18) 

Expanding the above formula about the fixed point c*, we have the following scaling matrix M (c*): 

BV 

BV 

M an (c *) = /-i(c*)-^(c*) 
09 

M na ( C *) = ^g( c *) 

M„„(c*)= (e+^(c*)J , (19) 

where a and j3 run from 1 to n — 1. Since V a (g) is a component perpendicular to g, one finds that V a (kc*) = for all 
k with help of Eqs.© and (USJ). This means that B g V a (c*) = 0. On the other hand, B g V n (c*) = 2V n (c*) /g* = -2e 
because V n (g) is quadratic in g. Therefore, 

M Q „ = 0, M„„ = -e (20) 
in Eq. I|19|) . Furthermore, we can rewrite M a( 3 (c*) in terms of /i Q ^ (a*). In fact, (a*) in Eq. Ijllfl is written as 

/ *x r 1 / *x 1 / *N 



Employing the following scaling properties: 



C(kg) = C(g) 

f a (kg) = kf a (g) 

dV a 2 dV a 

(kg) = k—(g), (22) 
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we get 

M a p (c*) = efi a (a*) . 
Eas. l(2"U|) and shows that M (c*) has a form of 



(23) 



M (c*) 








e/x (a*) 









* • ■ ■ * 


-e 



(24) 
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in the polar coordinates. It readily follows from this formula that M(c*)e n — — ee„. Thus we can derive all the 
eigenvalues of fJ.(a*) from M(c*l by removing — e, which is the eigenvalue corresponding to the eigenvector e„, from 
the set of the eigenvalues of M(c*), and, by multiplying by 1/e, the remaining eigenvalues. Further, if all the 
eigenvalues of M (c*) are negative, g(t, ao) behaves as 



9(t,a ) 



C(a*)t et 

according to Eq. (|13|l and the scaling property of C (a*) in Eq. 1221) . 



(25) 



IV. EXAMPLE 



Here are several examples. The first example is taken from the two-dimensional classical XY modelQ. Here, the 
beta function V (g) is given as 

vw =(_~i). <» 

for gi, (72 > 0. The deformed RGE has the fixed point c* = (e, e). The scaling matrix M (c*) of the deformed RGE is 
easily computed in terms of the cartesian coordinates as 

M(c*)= (_ e £ " 2£ ). (27) 

It has the eigenvalues — e and 2e. Employing Eq. I|12f> . we get 

which is a well-known result. As we have explained in the previous section, the other eigenvalue, — e, always appears 
in a deformed RGE © , which corresponds to the eigenvector c* /c* . 

The next example is the RGE in a one-dimensional quantum spin chain, studied by Itoi and Katop^|: it is defined 

by 

V(9)=(^]t 2 "±). (29) 



- ff2 (2 5l +A 52 ) 

The deformed RGE has the following three nontrivial fixed points: 



4,0),c-(0,±),c-^(-l,l). (30) 



The corresponding scaling matrices are 



*-l o M e ).^ = (_V-° e 

M 3 = ( V[ . (31) 

2-N 2-N J 
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The eigenvalues of those matrices are, respectively, 

N + 2 N + 2 2 + N 

—e, —c, and —e, (32) 

up to the common eigenvalue — e. The other eigenvalues divided by e are equal to those of the scaling matrices derived 
from the new RGE Q, which is computed in Ref.yf- It should be noted that the deformed RGEs in the above two 
examples do not correspond to those in 2 — e and 1 — e dimensions respectively. However, the derivation presented 
here is much simpler than the method using Eq.(J7J. 

The last example is the RGE in a field-theoretical model for nematic elastomers, proposed in Ref.0. In contrast 
to the previous examples, the deformed RGE is obtained exactly in 3 — e dimensions with 



v (a) = - l ( 9i (40.g? + 68. 9l .g 2 + 13«? 2 2 ) 
Ky > 8 (4 fl i + g 2 ) \ 2 52 (4«? 2 + 3 2 5l52 + 7.g 2 ) 



(33) 



Although V (g) is not quadratic polynomial, our result is applicable because all we need to apply the present method 
is the scaling property of V (g), Eq. ©. The deformed RGE has the three fixed points 

4-(±A <*) 



, 5 J " \59 59 y J V 7 . 

One can check that the scaling matrices have the following respective eigenvalues 

4e/5, -4e/59, and e/14 (35) 

in addition to the common eigenvalue — e. Now we turn to the case of just three dimensions. If <7i,<? 2 > 0, infrared 
behavior of a system is governed by the fixed point c^jSJ. Since the eigenvalue at is negative, g (t, a ) behaves as 

1 1 / 4 32\ 
ff ( t>ao )~-c a = T ^- > -j (36) 

for sufficiently large t, according to Eq. (|25p. The result is consistent with that in Ref.0. 



V. SUMMARY 



We have shown how to derive asymptotic behavior of a solution of RGE for infinite-order phase transition, by 
adding a linear term to this RGE. This method can allow us to apply a result of the e expansion to the case where 
e = Q. 
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